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Abstract 

We study a heavy piston that separates finitely many ideal gas particles moving 
inside a one-dimensional gas chamber. Using averaging techniques, we prove precise 
rates of convergence of the actual motions of the piston to its averaged behavior. The 
convergence is uniform over all initial conditions in a compact set. The results extend 
earlier work by Sinai and Neishtadt, who determined that the averaged behavior is 
periodic oscillation. In addition, we investigate the piston system when the particle 
interactions have been smoothed. The convergence to the averaged behavior again 
takes place uniformly, both over initial conditions and over the amount of smoothing. 

Mathematics Subject Classification (2000): 34C29, 37A60, 82C22 
Keywords: adiabatic piston, averaging, ideal gas 

1 Introduction 

Consider the following simple model of a piston separating two gas chambers: A piston of 
mass M 3> 1 divides a cylindrical chamber in M. 3 into two halves. The piston is parallel to 
the two ends of the chamber and can only move in the normal direction. On either side of 
the piston there are a finite number of gas particles of unit mass. All of the gas particles are 
point particles that interact with the walls of the chamber and with the piston via elastic 
collisions. The interactions of the gas particles with the piston and the ends of the chamber 
are completely specified by their motions along the normal axis of the chamber. Thus, this 
system projects onto a system inside the unit interval consisting of a massive point particle, 
the piston, which interacts with the gas particles on either side of it. These gas particles 
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make elastic collisions with the walls at the ends of the chamber and with the piston, but 
they do not interact with each other: They form an ideal gas. 

This simple model is useful for investigating the adiabatic piston problem This well- 
known problem from physics concerns an adiabatic piston, i.e. one with no internal degrees of 
freedom, that separates two gas chambers. Initially, the piston is fixed in place, and the gas 
in each chamber is in a separate thermal equilibrium. At some time, the piston is no longer 
externally constrained and is free to move. One hopes for an ergodic theorem to show that 
eventually the system will come to a full thermal equilibrium, where each gas has the same 
pressure and temperature. If one existed, then the final pressure and temperature could 
be predicted using thermodynamics. However, whether the system will evolve to thermal 
equilibrium and the interim behavior of the piston are mechanical problems that are not yet 
explicitly resolved [Zj. We do not give a full investigation of these issues here, but we wish to 
emphasize that they cannot be addressed by thermodynamics alone. The interested reader 
should see some papers [HJEI] and the references therein for details about recent progress in 
this area. 

Gruber et al. [HI El EH have extensively studied the simple model above using the Boltz- 
mann equation, the Liouville equation, and numerical simulations. Making various assump- 
tions, they observed that the system evolves in at least two stages. First, the system relaxes 
deterministically and adiabatically toward mechanical equilibrium, where the pressures on 
either side of the piston are equal. In the second, much longer, stage, the piston drifts 
stochastically in the direction of the hotter gas, and the temperatures of the gases equili- 
brate. Because of the piston's stochastic fluctuations, heat is allowed to flow between the 
gases, and in this sense the piston is no longer adiabatic. It should be emphasized that the 
behavior of the piston depends strongly on the ratios rii/M, where n\ and n 2 are the number 
of gas particles on either side of the piston. In the thermodynamic limit M, ni, n 2 — > oo while 
ni/M and n 2 / 'M are held fixed, Gruber et al. concluded that in the first stage above, the 
piston performs damped oscillatory motion, where the damping is strong if ni/M, n 2 /M > 1, 
and weak if rii/M,n 2 /M < 1. Using kinetic theory, Crosignani et al. [5j had already derived 
similar equations describing this damped oscillatory motion for the adiabatic piston. 

Sinai J7j also investigated the simple model, but with averaging techniques that examine 
the limit where M — > 00 while the total energy of the system is bounded and ni and n 2 are 
fixed. He determined that the averaged behavior of the piston is periodic oscillation, with the 
piston moving inside an effective potential well whose shape depends on the initial position 
of the piston and the initial energies of the gases. Neishtadt and Sinai [TH] pointed out that 
a classical averaging theorem due to Anosov, proved for smooth systems, can be extended 
to this case. Their insight allows us to conclude that if we examine the actual motions 
of the piston with respect to the slow time r = t/M 1 / 2 as M — ► 00, then in probability 
(with respect to Riemannian volume) most initial conditions give rise to orbits whose actual 
motion is accurately described by the averaged behavior for r e [0, 1]. 

This paper proves that the actual motions do not deviate by more than 0(M~ 1/2 ) from 
the averaged behavior for r G [0,1], i.e. for t E [0, M 1 ' 2 ]. Furthermore, the size of the 
deviations is bounded, independent of the initial conditions. 

We also investigate the behavior of the system when the interactions of the gas particles 
with the walls and the piston have been smoothed, so that Anosov's theorem applies directly. 
Let 5 > be a parameter of smoothing, so that 5 = corresponds to the hard core setting 
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Figure 1: The piston system with n\ = 3 and n 2 = 4. Note that the gas particles do not 
interact with each other, but only with the piston and the walls. 

above. Then the averaged behavior of the piston is still a periodic oscillation, which depends 
smoothly on 5. We show that the deviations of the actual motions of the piston from the 
averaged behavior are again not more than 0(M~ 1 ^ 2 ) on the time scale M 1//2 . The size of 
the deviations is bounded uniformly, both over initial conditions and over the amount of 
smoothing. 



2 Statement of results 

2.1 The hard core piston problem 

Consider the system of n\ + n 2 + 1 point particles moving inside the unit interval indicated 
in Figure ^ One distinguished particle, the piston, has position X and mass M. To the 
left of the piston there are ri\ > particles with positions x\j and masses mi t j, 1 < j < ni, 
and to the right there are n 2 > particles with positions x 2 j and masses m 2 j, 1 < j < n 2 . 
These gas particles do not interact with each other, but they interact with the piston and 
with walls located at the end points of the unit interval via elastic collisions. We denote the 
velocities by dX/dt = V and dxij/dt = Vij. There is a method for transforming this system 
into a billiard system consisting of a point particle moving inside an (ni +n 2 + l)-dimensional 
polytope jl], but we will not use this in what follows. 

We are interested in the dynamics of this system when the numbers and masses of the 
gas particles are fixed, the total energy is bounded, and the mass of the piston tends to 
infinity. When M = oo, the piston remains at rest, and each gas particle performs periodic 
motion. More interesting are the motions of the system when M is very large but finite. 
Because the total energy of the system is bounded, MV 2 /2 < cnst, and so V = C(M -1 / 2 ). 
Set e = M" 1 / 2 , and let W = V/e, so that dX/dt = eW with W = 0(1). 

When e = 0, the system has rii+n 2 + 2 independent first integrals (conserved quantities), 
which we take to be X, W, and sy = \vij\, the speeds of the gas particles. We refer to 
these variables as the slow variables because they should change slowly with time when e 
is small, and we denote them by h = (X, W, Sij, s 2 j) G M ni+n2 + 2 . Let h £ (t,z) = h e (t) 
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denote the dynamics of these variables in time for a fixed value of e, where z represents the 
dependence on the initial condition in phase space. We usually suppress the initial condition 
in our notation. Think of h e (-) as a random variable which, given an initial condition 
in the 2{n\ + n 2 + l)-dimensional phase space, produces a piecewise continuous path in 
]R"i+"2+ 2 . These paths are the projection of the actual motions in our phase space onto a 
lower dimensional space. The goal of averaging is to find a vector field on M. ni+n2+2 whose 
orbits approximate h £ (t). 

The averaged equation 

Sinai [T7] derived 



d 
dr 



X 
W 



H(h) :- 



W 

j = l m l,3 s l,j _ 22j = l m 2,3 s 2,j 



X 



»ijW 
1-X 
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as the averaged equation (with respect to the slow time r = et) for the slow variables. We 
provide a heuristic derivation in Section 13.21 Sinai solved this equation as follows: From 
d\n(sx :j )/dr = -d\n(X)/dr, s ld (r) = s 1}j (0)X(0)/X(r). Similarly, s 2J (r) = s 2J (0)(l - 
X(0))/(1 -X(r)). Hence 

d 2 x ET=i mi^i.(o) 2 x(o) 2 ET=i ^'^(o) 2 (i - *(o)) 2 



dr 2 



X 3 



[i-xy 



and so (X, W) behave as if they were the coordinates of a Hamiltonian system describing a 
particle undergoing periodic motion inside a potential well. If we let E { = J2™Li m i,j s l,j/2 
be the kinetic energy of the gas particles on one side of the piston, the effective Hamiltonian 
may be expressed as 



1 2 ^(0)X(0) 2 

2 X 2 



£ 2 (0)(1-X(0))^ 

(i-xy 



(2) 



Hence, the solutions to the averaged equation are periodic for all initial conditions under 
consideration. 

Main result in the hard core setting 

The solutions of the averaged equation approximate the motions of the slow variables, h e (t), 
on a time scale 0(1/ e) as e — > 0. Precisely, let h(r, z) = h(r) be the solution of 



dh 



H(h), h{0) = h e {0). 



Again, think of h(-) as being a random variable that takes an initial condition in our phase 
space and produces a path in R n i+ n 2+ 2 . 
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Next, fix a compact set V C R n ^+ n ^+ 2 such that h G V => X CC (0, 1), W CC R, and 
Sjj CC (0,oo) for each % and j. 1 For the remainder of this discussion we will restrict our 
attention to the dynamics of the system while the slow variables remain in the set V. To 
this end, we define the stopping time T e (z) = T e = inf{r > : h(r) ^ V or h e {r fe) V}. 

Theorem 2.1. For each T > 0, 

sup sup \h e (r/e) - h(r)\ = O(e) as e = M~ 1/2 -> 0. 

initial conditions 0<t<T/\T e 

s.t. h s (o)ev 

Note that the stopping time does not unduly restrict the result. Given any c such that 
h = c => X G (0, 1), Si j G (0, 00), then by an appropriate choice of the compact set V we 
may ensure that, for all e sufficiently small and all initial conditions in our phase space with 
h e (0) = c, T £ > T. We do this by choosing V 3 c such that the distance between dV and the 
periodic orbit h(r) with h(0) = c is positive. Call this distance al. Then T e can only occur 
before T if h e (r/e) has deviated by at least d from h(r) for some r G [0, T). Since the size 
of the deviations tends to zero uniformly with e, this is impossible for all small e. 

2.2 The soft core piston problem 

In this section, we consider the same system of one piston and gas particles inside the unit 
interval considered in Section CTT1 but now the interactions of the gas particles with the walls 
and with the piston are smooth. Let k : R — >• R be a C 2 function satisfying 

• k(x) — if x > 1, 

• k'(x) < if x < 1. 

Let 5 > be a parameter of smoothing, and set Kg(x) = k(x/8). Then consider the Hamil- 
tonian system obtained by having the gas particles interact with the piston and the walls 
via the potential 

^ Ks{xij) +Kg(X - Xij) + Ks{x 2 ,j ~ X)+K S (1- X 2 ,j). 
i=l 3=1 

As before, we set e = M' 1 / 2 and W = V/e. If we let 

E iJ = \ m ^-j v i j + K s( x i,j) + K s(X - xxj), 1 < j < Tlx, 

I ^ (3) 

E 2 ,j = 2 m 2,i^2,i + K s{x2,j ~ X) + K S (1 - X2j), 1 < j < n 2 , 

then Eij may be thought of as the energy associated with a gas particle, and W 2 /2 + 
Y^jLx E x,j + Sj=i E 2 j is the conserved energy. 

1 We have introduced this notation for convenience. For example, /i£ V=^X CC (0, 1) means that there 
exists a compact set A c (0, 1) such that h E V ^ X E A, and similarly for the other variables. 
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When e = 0, the Hamiltonian system admits rii + n 2 + 2 independent first integrals, 
which we choose this time as h = (X, W, E 1 j, E 2 j). While discussing the soft core dynamics 
we use the energies Eij rather than the variables Sjj = ^j2Ei ) j/m i ^, which we used for the 
hard core dynamics, for convenience. 

For comparison with the hard core results, we formally consider the dynamics described 
by setting 5 = to be the hard core dynamics described in Section 12.11 This is reasonable 
because we will only consider gas particle energies below the barrier height k(0). Then for 
any e, 5 > 0, h 5 £ (t) denotes the actual time evolution of the slow variables. While discussing 
the soft core dynamics we often use 5 as a superscript to specify the dynamics for a certain 
value of 5. We usually suppress the dependence on 5, unless it is needed for clarity. 

Main result in the soft core setting 

We have already seen that when 5 = 0, there is an appropriate averaged vector field H° 
whose solutions approximate the actual motions of the slow variables, h £ (t). We will show 
that when 5 > 0, there is also an appropriate averaged vector field H 5 whose solutions still 
approximate the actual motions of the slow variables, h £ (t). We delay the derivation of H s 
until Section EIU 

Fix a compact set V C W^+ n ^+ 2 such that h E V => X CC (0,1), W CC R, and 
Eij CC (0, k(0)) for each % and j. For each e, 5 > we define the functions h s (-) and T £ on 
our phase space by letting h 5 {r) be the solution of 

— = H\h s ), h'{0) = h s £ (0), (4) 
and T £ 5 = inf{r > : h 5 (r) <£ V or h 6 e (r/e) <£ V}. 

Theorem 2.2. There exists Sq > such that the averaged vector field H 5 {h) is C 1 on the 
domain {(5, h) : < 5 < 5q, h G V}. Furthermore, for each T > 0, 

sup sup sup \h 5 £ {r/e) - h 5 (r)\ = 0{e) as e = M~ 1/2 -> 0. 

0<<5<<5o initial conditions 0<t<TAT s 
s.t. h*(0)£V 

As in Section |2~TI for any fixed c there exists a suitable choice of the compact set V such 
that for all sufficiently small e and 5, T £ > T whenever h £ (0) = c. 

2.3 Applications and generalizations 

Relationship between the hard core and the soft core piston 

It is not a priori clear that we can compare the motions of the slow variables on the time 
scale 1/e for 5 > versus 5 = 0, i.e. compare the motions of the soft core piston with the 
motions of the hard core piston on a relatively long time scale. It is impossible to compare 
the motions of the fast-moving gas particles on this time scale as e — > 0. As we see in Section 
the frequency with which a gas particle hits the piston changes by an amount 0(5) when 
we smooth the interaction. Thus, on the time scale 1/e, the number of collisions is altered 
by roughly 0(5/e), and this number diverges if 5 is held fixed while e — ► 0. 
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Similarly, one might expect that it is impossible to compare the motions of the soft 
and hard core pistons as e — > without letting 5 — > with e. However, from Gronwall's 
Inequality it follows that if h s (0) = h°(0), then sup 0<r<TAT * AT o \ h 5 (r) — h°(r)\ = 0(6). From 
the triangle inequality and Theorems 12.11 and 12.21 we obtain the following corollary, which 
allows us to compare the motions of the hard core and the soft core piston. 

Corollary 2.3. As e = M" 1 / 2 , 5-^0, 

sup sup sup \h* e (t)-h°(t)\=0(e)+0(S). 

ceV initial conditions 0<t<(TAT s AT°)/e 
s.t. h|(0)=c=A°(0) ~ E 

This shows that, provided the slow variables have the same initial conditions, 

sup \h 5 £ {t)-h £ {t)\=O{e) + O{5). 

0<t<l/e 

Thus the motions of the slow variables converge on the time scale 1/e as e, 5 — > 0, and it is 
immaterial in which order we let these parameters tend to zero. 

The adiabatic piston problem 

We comment on what Theorem 12.11 says about the adiabatic piston problem. The initial 
conditions of the adiabatic piston problem require that W(0) = 0. Although our system is 
so simple that a proper thermodynamical pressure is not defined, we can define the pressure 
of a gas to be the average force received from the gas particles by the piston when it is held 
fixed, i.e. P 1 = 2m uSi/H = 2E i/ x and P 2 = 2E 2 /(1 - X). Then if P^O) > P 2 (0), 

the initial condition for our averaged equation ([T} has the motion of the piston starting at 
the left turning point of a periodic orbit determined by the effective potential well. Up to 
errors not much bigger than M~ 1//2 , we see the piston oscillate periodically on the time scale 
M 1/2 . If Pi(0) < P 2 (0), the motion of the piston starts at a right turning point. However, if 
Pl(0) = p2(0), then the motion of the piston starts at the bottom of the effective potential 
well. In this case of mechanical equilibrium, h(r) = h(0), and we conclude that, up to errors 
not much bigger than M™ 1 / 2 , we see no motion of the piston on the time scale M 1//2 . A 
much longer time scale is required to see if the temperatures equilibrate. 

Generalizations 

There are several other ways this work could be generalized. For example, one could replace 
the walls at and 1 by heat baths, so that whenever a gas particle reaches a heat bath, 
it is absorbed, and another gas particle is emitted with its velocity chosen independently 
and randomly according to some distribution. Similarly, the walls could be replaced by a 
forcing mechanism (much like the bumpers in a pinball machine) that changes a colliding 
gas particle's kinetic energy. 

A simple generalization of Theorem I2.1[ proved by similar techniques, follows. The 
system consists of iV — 1 pistons, that is, heavy point particles, located inside the unit 
interval at positions X\ < X 2 < . . . < X^-i- Walls are located at X = and X^ = 1, 
and the piston at position X{ has mass M: L . Then the pistons divide the unit interval into iV 
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chambers. Inside the i chamber, there are rij > 1 gas particles whose locations and masses 
will be denoted by Xij and m it j, respectively, where 1 < j < rij. All of the particles are 
point particles, and the gas particles interact with the pistons and with the walls via elastic 
collisions. However, the gas particles do not directly interact with each other. We scale the 
piston masses as Mj = Mi/e 2 with Mi constant, define Wi by dXi/dt = eWi, and let Ei be 
the kinetic energy of the gas particles in the i th chamber. Then we can find an appropriate 
averaged equation whose solutions have the pistons moving like an N- dimensional particle 
inside a potential well with an effective Hamiltonian 

If we write the slow variables as h — (Xi, Wi, \vij\) and fix a compact set V such that 
h G V X i+ i — Xi CC (0,1), Wi CC R, and \vij\ CC (0, oo), then the convergence of 
the actual motions of the slow variables to the averaged solutions is exactly the same as the 
convergence given in Theorem 12.11 



3 Preparatory material 

Before the proofs of our main results, we present some averaging results, as well as a heuris- 
tic derivation of the averaged equation for the hard core piston. This material provides 
background for our work and establishes some notation. 

3.1 The averaging framework 

In this section, consider a family of ordinary differential equations 

| = Z(»,e) (5) 

on a smooth, finite-dimensional Riemannian manifold Ai, which depends on the real param- 
eter £ G [0, £o\. Assume 

• Regularity: the functions Z and dZ/de are both C 1 on M. x [0, e ]- 

We denote the flow generated by Z(-,e) by z e (t,z) = z £ {t). We will usually suppress the 
dependence on the initial condition z = z £ (0,z). Think of z E (-) as being a random variable 
whose domain is the space of initial conditions for the differential equation J5J) and whose 
range is the space of continuous paths (depending on the parameter t) in M.. 

• Existence of smooth integrals: zo(t) has m independent C 2 first integrals h = (hi, . . . , h m ). 

Then h is conserved by Zo(t), and at every point the linear operator dh/dz has full rank. It 
follows from the implicit function theorem that each level set M. c = {h = c} is a smooth 
submanifold of co-dimension m, which is invariant under zo(t). Further, assume that there 
exists an open ball U C M. m satisfying: 
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eH(z £ ,e) 



• Compactness: Vc G U, A4 C is compact. 

• Preservation of smooth measures: Vc G U, zq(1;)\m c preserves a smooth measure fi c 
that varies smoothly with c, i.e. there exists a C 1 function / : M. — > IR>o such that 
J\m c is the density of /i c with respect to the restriction of Riemannian volume. 

Set h e (t) = h(z £ (t)). Since dh /dt = 0, Hadamard's Lemma allows us to write 

dh £ 
~dt 

for some function H : M. x [0, £q] W. Observe that 
dh 

-^-(t) = Dh(z £ (t))Z(z £ (t) } e) = Dh(z £ (t))(Z(z £ (t),e) - Z(z £ (t),0)), 

so that 

H(z,0)=£sz ]e=o h. 

Here £ denotes the Lie derivative. 

Define the averaged vector field H by 



H(h)= / H(z,0)dfx h (z). (6) 

JMh 

Then H is C 1 . Fix a compact set VcW, and let h{r) be the solution of 

dh - - 

- = H(h), h(0) = h £ (0). 
dr 

We only consider the dynamics in a compact subset of phase space, so for initial conditions 
z G h~ l U, define the stopping time T £ = inf{r > : h(r) ^ V or h £ {r/e) ^ V}. 

Heuristically, think of the phase space M. as being a fiber bundle whose base is the 
open set U and whose fibers are the compact sets M.h- See Figure 121 Then the vector 
field Z(-,0) is perpendicular to the base, so its orbits z (t) flow only along the fibers. Now 
when 0<£<1, the vector field Z(-,e) acquires a component of size 0(e) along the base, 
and so its orbits z £ (t) have a small drift along the base, which we can follow by observing 
the evolution of h £ (t). Because of this, we refer to h as consisting of the slow variables. 
Other variables, used to complete h to a parameterization of (a piece of) phase space, are 
called fast variables. Note that h £ (t) depends on all the dimensions of phase space, and so 
it is not the flow of a vector field on the m-dimensional space U. However, because the 
motion along each fiber is relatively fast compared to the motion across fibers, we hope to 
be able to average over the fast motions and obtain a vector field on IA that gives a good 
description of h £ {t) over a relatively long time interval, independent of where the solution 
z £ (t) started on M.h e (o)- Because our averaged vector field, as defined by Equation (jUJ), only 
accounts for deviations of size 0(e), we cannot expect this time interval to be longer than 
size 0(l/e). In terms of the slow time r = et, this length becomes 0(1). In other words, 
the goal of the first-order averaging method described above should be to show that, in some 
sense, sup 0<T<1ATe \h £ (r/e) — h(r)\ — ► as e — > 0. This is often referred to as the averaging 
principle. 
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Figure 2: A schematic of the phase space Ai. 
Some averaging results 

So far, we are in a general averaging setting. Frequently, one also assumes that the invariant 
submanifolds, A4h, are tori, and that there exists a choice of coordinates z = (h, ip) on Ai 
in which the differential equation © takes the form 

^=eH{h,ip,£), ( ^- = ${h,ip,e). 

Then if ip e S 1 and the differential equation for the fast variable is regular, i.e. ip, 0) is 
bounded away from zero for h G U, 

sup sup \h e (r/e) — h(r) \ = 0{e) as e — > 0. 

initial conditions 0<t<1AT e 

s.t. h e (o)eV 

See for example Chapter 5 in ^1 or Chapter 3 in |13j . 

When the differential equation for the fast variable is not regular, or when there is more 
than one fast variable, the typical averaging result becomes much weaker than the uniform 
convergence above. For example, consider the case when ip e T n , n > 1, and the unperturbed 
motion is quasi-periodic, i.e. <&(h,<p,0) = Q(h). Also assume that H £ C n+2 and that Q is 
nonvanishing and satisfies a nondegeneracy condition on U (for example, Q : IA — > T n is a 
submersion). Let P denote Riemannian volume on .M. Neishtadt fT3*l ITij showed that in 
this situation, for each fixed 5 > 0, 

P ( sup \h E (r/e) - h(r)\ > 6) = 0(y/e/6), 

\0<t<1AT £ J 

and that this result is optimal. Thus, the averaged equation only describes the actual motions 
of the slow variables in probability on the time scale 1/e as e — > 0. 
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Neishtadt's result was motivated by a general averaging theorem for smooth systems due 
to Anosov. This theorem requires none of the additional assumptions in the averaging results 
above. Under the conditions of regularity, existence of smooth integrals, compactness, and 
preservation of smooth measures, as well as 

• Ergodicity: for Lebesgue almost every c G U, (zo(-),fi c ) is ergodic, 

Anosov jl| HSj showed that 

sup \h £ (r/e) — h(r)\ — > in probability (w.r.t. Riemannian volume) as e — > 0. 

0<t<1AT £ 

If we consider h e (-) and h(-) to be random variables, this is a version of the weak law 
of large numbers. In general, we can do no better: There is no general strong law in this 
setting. There exists an example due to Neishtadt (which comes from the equations for the 
motion of a pendulum with linear drag being driven by a constant torque) where for no 
initial condition in a positive measure set do we have convergence of h e (t) to h(et) on the 
time scale 1/e as e — > Here, the phase space is M X S 1 , and the unperturbed motion 
is (uniquely) ergodic on all but one fiber. 

In 2004, Kifer ^T] gave necessary and sufficient conditions for the averaging principle to 
hold in an averaged sense with respect to initial conditions. He also showed explicitly that 
his conditions are met in the setting of Anosov's theorem. A recent, relatively simple proof 
of the version of Anosov's theorem stated above is due to Dolgopyat jO]. 



3.2 Heuristic derivation of the averaged equation for the hard core 
piston 

We present here a heuristic derivation of Sinai's averaged equation (JIJ that is found in jH]. 

First, we examine interparticle collisions when e > 0. When a particle on the left, say the 
one at position xij, collides with the piston, sij and W instantaneously change according 
to the laws of elastic collisions: 



V+ 



m lj + M 



mi j - M 



2m 



1,3 



2M 
M -m 1 









y- 



(7) 



If the speed of the left gas particle is bounded away from zero, and W = M X I 2 V is also 
bounded, it follows that for all e sufficiently small, any collision will have v^j > and 
< 0. In this case, when we translate Equation (J2J) into our new coordinates, we find that 



%3 



so that 



W+ 



1 



1 + e mij 



1 — s 2 mij 



-2e 
1 — e 2 m\ 









w- 



(8) 



AW 



S l,3 S l,3 

w + - w~ 



-2eW~ + 0(e 2 ), 
= +2emijs^. + 0(e 2 
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The situation is analogous when particles on the right collide with the piston. For all e 
sufficiently small, s 2 j and W instantaneously change by 

AW = W + -W- = -2em 2jj s 2j + 0(e 2 ), 

As 2,j = 4,j ~ «2,j = +2eW~ + 0(e 2 ). 

We defer discussing the rare events in which multiple gas particles collide with the piston 
simultaneously, although we will see that they can be handled appropriately 

Let At be a length of time long enough such that the piston experiences many collisions 
with the gas particles, but short enough such that the slow variables change very little, 
in this time interval. From each collision with the particle at position xxj, W changes by 
an amount +2emx t jSxj + 0(e 2 ), and the frequency of these collisions is approximately 
Arguing similarly for collisions with the other particles, we guess that 

l\W v-^ Si 7 - v-^ S 2 7 9n 

— = e £ Jm^ - e £ 2^,^^ + 

3=1 3=1 V ' 

Note that not only does the position of the piston change slowly in time, but its velocity 
also changes slowly, i.e. the piston has inertia. With r = et as the slow time, a reasonable 
guess for the averaged equation for W is 

dW _ ET=i m U s lj _ E?=i m 2,j4,j 
dr X l-X 

Similar arguments for the other slow variables lead to the averaged equation ((TJ). 



4 Proof of the main result for the hard core piston 

4.1 Proof of Theorem 12.11 with only one gas particle on each side 

We specialize to the case when there is only one gas particle on either side of the piston, 
i.e. we assume that n\ = n 2 = 1. We then denote x^i by X\, m 2i2 by 7712, etc. This allows the 
proof's major ideas to be clearly expressed, without substantially limiting their applicability. 
At the end of this section, we outline the simple generalizations needed to make the proof 
apply in the general case. 



A choice of coordinates on the phase space for a three particle system 

As part of our proof, we choose a set of coordinates on our six-dimensional phase space such 
that, in these coordinates, the e = dynamics are smooth. Complete the slow variables 
h = (X, W, Si, s 2 ) to a full set of coordinates by adding the coordinates ^ G [0, 1]/ ~ 1 = 
S 1 , i = 1, 2, defined as follows: 



<pi = ipiix^v^X) 



2X 



if v x > 
1 - f£ if v x < 



I 9TT^Y if^2<0' 

1 " W^t) ^ ^2 > 
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When e = 0, these coordinates are simply the angle variable portion of action-angle co- 
ordinates for an integrable Hamiltonian system. They are defined such that collisions 
occur between the piston and the gas particles precisely when ipi or p 2 — 1/2. Then 
z = (h, (pi, (p 2 ) represents a choice of coordinates on our phase space, which is homeomorphic 
to (a subset of M 4 ) x T 2 . We abuse notation and also let h(z) represent the projection onto 
the first four coordinates of z. 

Now we describe the dynamics of our system in these coordinates. When (pi, p 2 7^ 1/2, 



d<pi 
dt 

dp 2 
~dT 



2X X ^1 
gl I eW 
2X X 



S2 



2(1-X) 

S2 



if < ip x < 1/2 
\-<p\) if 1/2 < y?i < 1 

+ f^p 2 if < <^ 2 < 1/2' 

-^(1-^) if 1/2 < ^ 2 < 1 



2(1-X) 

Hence between interparticle collisions, the dynamics are smooth and are described by 



eW, 



dX 
~dt 

dW dsi ds 2 
dt dt dt 

d<fi Si 



dt 

dip 2 

~dT 



2k + °^ 



(9) 



«2 



2(1 -X) 



+ 0(e) 



When pi reaches 1/2, while p 2 / 1/2, the coordinates X, s 2 ,ifi, and (p 2 are instanta- 
neously unchanged, while Si and W instantaneously jump, as described by Equation (jHJ). It is 
curious that sf +eW + = s± —eW~, so that dipi/dt is continuous as (fx crosses 1/2. However, 
the collision induces discontinuous jumps of size 0(e 2 ) in dX/dt and dp 2 /dt. Denote the 

1 -2e 
2em\ 1 

The situation is analogous when p 2 reaches 1/2, while (pi / 1/2. Then W and s 2 are 

1 —2em 2 
2e 1 



linear transformation in Equation (jHJ) with j — 1 by A ltE . Then A lfE 



+ 0(e 2 



instantaneously transformed by a linear transformation A 2j£ 

We also account for the possibility of all three particles colliding simultaneously. There 
is no completely satisfactory way to do this, as the dynamics have an essential singularity 
near {pi = p 2 = 1/2}. Furthermore, such three particle collisions occur with probability 
zero with respect to the invariant measure discussed below. However, the two 3x3 matrices 



~A 1;£ 




"1 " 


1 




A 2i£ _ 



have a commutator of size 0(e 2 ). We will see that this small of an error will make no 
difference to us as e — > 0, and so when p x = <p 2 = 1/2, we pretend that the left particle 
collides with the piston instantaneously before the right particle does. Precisely, we transform 
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the variables si, W, and s 2 by 



We find that 



W 

st 



AW = W 

As 2 = s% 



1 " 




~A l)E 0" 


A 2 , £ _ 




1 



°1 



-2eW~ + 0{e 2 ), 

= +2emis^ — 2em 2 S2 + 0(e 2 ), 

+2eW- + 0(e 2 ). 



The above rules define a flow on the phase space, which we denote by z £ (t). We denote 
its components by X E (t), W e (t), si )E (t), etc. When e > 0, the flow is not continuous, and for 
definiteness we take z £ {t) to be left continuous in t. 

Because our system comes from a Hamiltonian system, it preserves Liouville measure. In 
our coordinates, this measure has a density proportional to X(l — X). That this measure 
is preserved also follows from the fact that the ordinary differential equation Q preserves 
this measure, and the matrices A% je , A 2fi have determinant 1. Also note that the set {ipi = 
ip 2 = 1/2} has co-dimension two, and so {J t z e (t){ipi = ip 2 = 1/2} has co-dimension one, 
which shows that only a measure zero set of initial conditions will give rise to three particle 
collisions. 



Argument for uniform convergence 

Step 1: Reduction using Gronwall's Inequality. Define H(z) by 

W 

2m 1 s 1 5 ipi=1/2 - 2m 2 s 2 S ip2=1 / 2 



Hz) 



-2W5, 



2W5, 



¥>i=l/2 



<P2=l/2 



Here we make use of Dirac delta functions. All integrals involving these delta functions 
may be replaced by sums. We explicitly deal with any ambiguities arising from collisions 
occurring at the limits of integration. 

Lemma 4.1. For < t < 

— — £ ' 

h e (t)-h e (0)=e f H(z £ (s))ds + 0(e), 
Jo 

where any ambiguity about changes due to collisions occurring precisely at times and t is 
absorbed in the 0(e) term. 

Proof. There are four components to verify. The first component requires that X £ (t) — 
X £ (0) = e jj W £ (s)ds + 0(e). This is trivially true because X £ {t) - X £ (0) = e f* W £ (s)ds. 



The second component states that 



W £ (t) - W £ (0) = e 2m 1 s 1>£ (s)5 (pi ,( s )=i/ 2 - 2m 2 s 2i£ (s)(5 (/P2 ^ s)=1 / 2 ds + 0(e). 
Jo 



(10) 
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Let r k and qj be the times in (0, t) such that y?i, e (r fc ) = 1/2 and (p 2 ,e{qj) — 1/2, respectively. 
Then 

W £ (t) - W e (0) = J2 &W £ (r k ) + £ AW e fe) + 0(e). 

Observe that there exists u; > such that for all sufficiently small e and all h G V, 1/u; < 
^ < uj. Thus the number of collisions in a time interval grows no faster than linearly in 
the length of that time interval. Because t < T/e, it follows that 

W £ (t) -W £ (0)=eJ2 2m lSh£ {r k ) - ej^ 2m 2 s 2 , 6 ( qj ) + 0(e), 

and Equation ()10|) is verified. Note that because V is compact, there is uniformity over all 
initial conditions in the size of the 0(e) terms above. The third and fourth components are 
handled similarly. □ 



Next, h(r) satisfies the integral equation 



while h £ {r/e) satisfies 



h(r) - h(0) 



h £ (r/s)-h £ (0) = O(e)+e [* H 

Jo 



H(h(a))da, 



z £ (s))ds 



r/e 



0(e) +e / H(z £ (s)) - H(h £ (s))ds + H(h £ (a/e))da 



for < r < T A T F 



Define e £ (r) = e f^ £ H(z £ (s)) — H(h £ (s))ds. It follows from Gronwall's Inequality that 



sup \h(r) - h £ (T/e)\ < 0(e) + sup |e e (r)| ) e 

0<r<TAT, V 0<t<TAT, 



Lip(H[ v )T 



(11) 



Gronwall's Inequality is usually stated for continuous paths, but the standard proof (found 
in jTB]) still works for paths that are merely integrable, and \h(r) — h £ {j/e)\ is piecewise 
smooth. 



Step 2: A splitting according to particles. Now 



2miSi^ 1= i/ 2 - misj/X 
-2W5 vl=1/2 + Sl W/X 




H(z) - H(h) 



+ 







-2m 2 S2^ 2 =i/2 + m 2 sl/ (1 - X) 


2W6< P2=1/2 -s 2 W/(l-X) 
and so, in order to show that sup 0<r<TATe |e e (r)| = 0(e), it suffices to show that 

Sl,e(s)S tpie ( s ) = i/ 2 



sup 

0<T<TAT e 



sup 

0<t<T/\T e 



2XJs) 



-ds 



^(l), 



^ £ (g)gl, £ (g) 



ds 



Oil), 
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as well as two analogous claims about terms involving ip 2 ,e- Thus we have effectively separated 
the effects of the different gas particles, so that we can deal with each particle separately. 



We will only show that sup 0<T<TATe 
three terms can be handled similarly. 



pi, e (s)=l/2 2X £ (s) 



ds 



0(1). The other 



Step 3: A sequence of times adapted for ergodization. Ergodization refers to the 
convergence along an orbit of a function's time average to its space average. Because of the 
splitting according to particles above, one can easily check that | f Q H(zo(s))ds = H(h ) + 
0(l/t), even when zq(-) restricted to the invariant tori M.h Q is not ergodic. In this step, for 
each initial condition z £ (0) in our phase space, we define a sequence of times t^^ inductively 
as follows: t 0£ = inf{t > : y?i,e(£) = 0}, t k+ ^ £ = inf{t > £ fc)£ : y2i, e (£) = 0}. This sequence is 
chosen because ( s )=i/2 is "ergodizd" as time passes from tk,o to tfc+i,o- If £ is sufficiently 
small and tk+i )E < (T A T £ )/e, then the spacings between these times are uniformly of order 
1, i.e. 1/u < t k +i,e - t kj£ < uj. Thus, 



sup 

Q<r<TAT e 



r/e 



Sl,e(s)5 lf>1 Js)=l/2 ~ 



2XJs) 



-ds 



<o(i)+ J2 



tk+l,e 



tk,e 



2X £ (s) 



-ds 



(12) 



Step 4: Control of individual terms by comparison with solutions along fibers. 

The sum in Equation (|T2*|) has no more than 0(1/ e) terms, and so it suffices to show that 
each term is no larger than 0(e). We can accomplish this by comparing the motions of z £ (t) 
for t k)£ <t< tfc + i )£ with the solution of the e = version of Equation Q that, at time tk, £ , is 
located at z £ (tk, £ )- Since each term in the sum has the same form, without loss of generality 
we will only examine the first term and suppose that io,e — 0, i.e. that (fi^ £ (0) = 0. 

Lemma 4.2. If t 1)£ < then sup < t < tle |*b(*) - z £ (t)\ = 0(e). 

Proof. To check that sup 0<t<tl£ \h (t) — h £ (t)\ = 0(e), first note that ho(t) = h (0) = h £ (0). 
Then dX £ /dt = 0(s), so that X (t) — X £ {t) = 0(et). Furthermore, the other slow variables 
change by 0(e) at collisions, while the number of collisions in the time interval [0, ti £ ] is 
Oil). 

It remains to show that sup 0<t<ti |y>i,o(0 — fi,e(t)\ — 0(e). Using what we know about 
the divergence of the slow variables, 

" = jf %$T) ~ + ° {F - )is = SI ° {F - )ia = ° (£) 

for < t < ti j£ . Showing that sup Q<t<t |v ? 2,o(^) ~ ^2,e(t)\ = 0(e) is similar. 

□ 
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From Lemma t 1>£ = t lfi + 0(e) = 2X /s lfi + 0(e). We conclude that 



si, £ (s)^ li£ ( s )=i/2 - 77Trj-rds = 0(e) + / s lfi (s)5^ lAs)=1/2 



2X £ {s) 



2X (s] 



ds 



0(e) + si, - ti, £ 



J i,o 

2x; 



0(e). 



From Equations (|TT|) and (JT2J), we see that sup 0<r<TAT£ \h £ (r/e) — h(r)\ = 0(e), inde- 
pendent of the initial condition in h~ x V. 



4.2 Extension to multiple gas particles 

When ni,n 2 > 1, only minor modifications are necessary to generalize the proof above. We 
start by extending the slow variables h to a full set of coordinates on phase space by defining 
the angle variables ipij G [0, 1]/ ~ 1 = S 1 for 1 < i < 2, 1 < j < nf. 



<f2,j = r->.j(-'\j- <2, r -V) 



1—X2- 



2(1- V ) 



1-x 



2.J 



2(1-X) 



if u 2)3 - < 
if v 2 ,j > 



Then d(pij/dt = s hj {2X)~ 1 +0(e), d(p 2J /dt = s 2 , j (2(l-X))- 1 + 0{e), and z = (h,(pij,(p 2 J) 
represents a choice of coordinates on our phase space, which is homeomorphic to (a subset of ]R n i+ n 2+ 2 ) x 
qpni+n 2 _ j n these coordinates, the dynamical system yields a discontinuous flow z £ (t) on phase 
space. The flow preserves Liouville measure, which in our coordinates has a density propor- 
tional to X Ul (l — X) TL2 . As is Section l4*7T) one can show that the measure of initial conditions 
leading to multiple particle collisions is zero. 
Next, define H(z) by 



H(z) 



W 



EJli 2m ld s ld 5^ =1/2 - Yj%i 2m 2 s 2J 5 V2> . =1/2 
-2Wc 
2WS t 



-2W5 91d=1/2 



><P2,j=l/2 



For < t < x ^ k , h £ {t) - h e (0) = e f*H(z £ (s))ds + 0{e). From here, the rest of the proof 
follows the same arguments made in Section 14.11 



5 Proof of the main result for the soft core piston 

For the remainder of this work, we consider the family of Hamiltonian systems introduced in 
Section|2I21 which are parameterized by e, 5 > 0. For simplicity, we specialize to n\ = n 2 = 1. 
As in Section HJ the generalization to rii,n 2 > 1 is not difficult. The Hamiltonian dynamics 
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are given by the following ordinary differential equation: 



dX 

~dt 
dW 

~dT 

dxi 

~dT 



eW, 

e (-k' s (X - Xl ) + k' 5 (x 2 - X)) 

Vl, 



dX2 

~dT 



(13) 



«2, 



^7 = — (-«J(x 2 -X) + 4(1 -x 2 )). 
dt m 2 

Recalling the particle energies defined by Equation (JHJ), we find that 

^ = eW^X - xi), ^ = -£H/<(x 2 - X). 

For the compact set V introduced in Section 12.21 fix a small positive number £ and 
an open set U C K 4 such that V C W and h E U ^ X E (£,1 — £), W CC R, and 
£ < Ei, E 2 < k(0) — £. We only consider the dynamics for < 5 < £/2 and h eW. 

Define 

Lq(x x ) = £/i(xi,X, <5) = k s (xi) + /^pf - xi), 
U 2 {x 2 ) = U 2 (x 2 , X, 5) = k & (x 2 - X) + Kg(l - x 2 ). 

Then the energies Ei satisfy Ei = rriivf /2 + Ui(xi). 

Let Ti = Ti(X, Ei, S) and T 2 = T 2 (X, E 2 , S) denote the periods of the motions of the left 
and right gas particles, respectively, when e = 0. 



Lemma 5.1. Fori = 1,2, 

Ti e C X {{X, E u 5): X e {£,!-£), E t e (£, k(0) - £ ), < 5 < £/2}. 
Furthermore, 



Ti{X,E u S) = J^X + 0(5), T 2 (X,E 2 ,6) = J^(l-X) + 0(5). 

The proof of this lemma is mostly computational, and so we delay it until Section El 
Note especially that the periods can be suitably defined such that their regularity extends 
to 5 = 0. 

In this section, and in Section |H1 below, we adopt the following convention on the use of 
the O notation. All use of the O notation will explicitly contain the dependence on e and 
5 as e, 5 — > 0. For example, if a function f{h,e,5) = 0(e), then there exists 5',e' > such 
that swpQ <£ < £ , !0< s<s',hev l/(M,<*)/e| < 00. 
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When e = 0, ^ = ±^£(£, - ^(a*))- Define a = a^, <5) by «*(a) = K{a/S) = E h 

so that a(E\,8) is a turning point for the left gas particle. Then a = 5n~ l {E i ) 1 where 
is defined as follows: k : [0, 1] — > [0, k(0)] takes to k(0) and 1 to 0. Furthermore, 
k G C 2 ([0, 1]), k' < 0, and k'(x) < if x < 1. By monotonicity, k _1 : [0, re(0)] — > [0, 1] exists 
and takes to 1 and k(0) to 0. Also, by the Implicit Function Theorem, k,^ 1 G C 2 ((0, k(0)]), 
< for y > 0, and — > — oo as y — »■ + . Because we only consider energies 
Ei G (£, k(0) — £), it follows that a{E i: 5) is a C 2 function for the domains of interest. 



5.1 Derivation of the averaged equation 

As we previously pointed out, for each fixed 5, Anosov's Theorem applies directly to the 
family of ordinary differential equations in Equation (jlHj) . provided that 5 is sufficiently 
small. The invariant fibers Aih of the e = flow are tori described by a fixed value of the 
four slow variables and {(X, W, x±, vi, X2, v 2 ) : E\ = mivf/2 + Ui(xi,X,5), E 2 = m 2 v 2 /2 + 
U2(x2, X, 5)}. If we use (xi,x 2 ) as local coordinates on M. h) which is valid except when 
v\ or t> 2 = 0, the invariant measure \ih of the unperturbed flow has the density 

dx\dx2 



T ly [±{Ey - tM*!)) T 2 ^ 2 (E 2 - U 2 {x 2 )) 



The restricted flow is ergodic for almost every h. See Corollary 16. II in Section 
Now 

W 

-k' s {X - xi) + k' s {x 2 - X) 
Wk' s (X - xt) 
-Wk' s {x 2 -X) 

and 



dj4 

dt 



k' s (X - x 1 )dyL h 



T-l 



X-a 



\ / 2m 1 



dx\ 



Ex 



k' s (X - x x ) 



Ti 



du 



y/Ei — U 



X-a 



dx\ 



k' s (X - Xi 



X-6 



Ei - k 5 (X - xi) 



Similarly, 



/ k' s (x 2 -X)dfi h 



^8m 2 E 2 
To ' 



It follows that the averaged vector field is 



H\h) 



W 



v^mi-Ei v / 8m2-E2 

Ti T 2 

TT/ \ / 8m 1 Ei 



T2 
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where from Lemma I3TT1 we see that H'(-) E C 1 ({(<5, h) : < S < 5/2, h E V}). H°(h) agrees 
with the averaged vector field for the hard core system from Equation (JTJ, once we account 
for the change of coordinates Ei = mjsf/2. 

Remark 5.1. An argument due to Neishtadt and Sinai ^H] shows that the solutions to the 
averaged equation (J3J) are periodic. This argument also shows that, as in the case 6 = 0, the 
limiting dynamics of (X, W) are effectively Hamiltonian, with the shape of the Hamiltonian 
depending on 5, X(0), and the initial energies of the gas particles. The argument depends 
heavily on the observation that the phase integrals 



I l (X,E l ,S) = 




dxdv 

<2+Ui(x,X,5)<Ei 



are adiabatic invariants, i.e. they are integrals of the solutions to the averaged equation. 
Thus the four-dimensional phase space of the averaged equation is foliated by invariant two- 
dimensional submanifolds, and one can think of the effective Hamiltonians for the piston as 
living on these submanifolds. 



5.2 Proof of Theorem 12^21 

The following arguments are motivated by our proof in Section HJ although the details are 
more involved as we show that the rate of convergence is independent of all small S. 



A choice of coordinates on phase space 



We wish to describe the dynamics in a coordinate system inspired by the one used in Section 
14.11 For each fixed S E (0, So], this change of coordinates will be C l in all variables on the 
domain of interest. However, it is an exercise in analysis to show this, and so we delay the 
proofs of the following two lemmas until Section |U1 



We introduce the angular coordinates ipi E [0, 1]/ 



(pi = p 1 (x 1 ,v 1 ,X) 



"XI 



ip 2 = ip 2 {x 2 ,v 2 ,X) 





j\ L 
1/2 

i-j-r 

Ti Ja 



- f 

T 2 Jx 2 

1/2 



mi/2 
Eh.-Ui(s) 



xi 



-1-a 



m 2 /2 
E 2 -U 2 (s) 



1-a 
T 2 Jx 2 



S 1 defined by 

if X\ = a 
ds if Vi > 

if xi = X — a 
ds if Vi < 

if x 2 = 1 — a 
ds if v 2 < 

if x 2 = X + a 
ds if v 2 > 



mi/2 
Ei-Uiis) 



m 2 /2 
E 2 -U 2 (s) 



(14) 



Then z = (h,ipi,(p 2 ) is a choice of coordinates on hr x U. As before, we will abuse notation 
and let h(z) denote the projection onto the first four coordinates of z. 

There is a fixed value of So in the statement of Theorem 12. 21 However, for the purposes of 
our proof, it will be convenient to progressively choose <5 smaller when needed. At the end 
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of the proof, we will have only shrunk 5q a finite number of times, and this final value will 
satisfies the requirements of the theorem. Our first requirement on 5 is that it is smaller 
than S/2. 

Lemma 5.2. If 5q > is sufficiently small, then for each 5 G (0, 8q] the ordinary differential 
equation (JT3J) in the coordinates z takes the form 



(15) 



where Z 5 G C\hT x U x [0, oo)). When z G hr x U, 



Z s (z,s) 



eW 

s(-k' s (X - xi(z)) + k' s (x 2 (z) 
eWk' s (X - xxiz)) 
-£Wk' s {x 2 {z)-X) 
± + O(e) 
i + 0{e) 



X)) 



(16) 



Recall that, by our conventions, the 0(e) terms in Equation ()16j) have a size that can 
be bounded independent of all 5 sufficiently small. Denote the flow determined by Z s (-,e) 
by z*(t), and its components by X s £ {t), W 5 £ {t), E( >£ (t), etc. Also, set h 5 £ (t) = h(z 5 £ (t)). From 
Equation (JTHJ), 



H d (z,e) 



ldh s £ 
e dt 



W 

-k' s (X - xx{z)) + k' s (x 2 (z) 
Wk' 5 {X-x x {z)) 
-Wk' s (x 2 (z)-X) 



(17) 



In particular, H s (z,e) = H s (z,0). 

Before proceeding, we need one final technical lemma. 

Lemma 5.3. If 5q > is chosen sufficiently small, there exists a constant K such that for 
allSe (0,5 ], k' 5 {\X - Xi{z)\) =0 unless <p t G [1/2 - K5, 1/2 + K5] . 



Argument for uniform convergence 

We start by proving the following lemma, which essentially says that an orbit z £ (t) only 
spends a fraction 0(5) of its time in a region of phase space where \H S (z s £ (t) , s)\ — \tj8(~s 
is of size (9(5 _1 ) 



H 5 (z 5 £ (t)M 



Lemma 5.4. For < T < T < 



TAT} 



f 

Jt> 



H s (z s £ (s),0)ds = 0(1 V (T - T)). 
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Proof. Without loss of generality, T = 0. From Lemmas 15.11 and 15.21 it follows that if we 
choose Sq, sufficiently small, then there exists u > such that for all sufficiently small e and 

all 5 e (0,5 ], h e V ^ 1/cu < ^ < to. Define the set B = [1/2 - K8, 1/2 + K5), where 

K comes from Lemma 15.31 Then we find a crude bound on j Q T \k' s (X £ (s) — Xi(z £ (s))) | ds 
using that 

M,e . j > 1/w if G £ 

IS 



dt 



<uo if ^ e g B c . 



This yields 



o 



|«4(X e '( a )-xi(«f( a )))|da< 



< 



cnst 
cnst 



T 



2Ku5 



ds 



-T + 2KuS = O(IVT). 



5 l 2 K^ + i^ 

Similarly, J T |4(x 2 (z£(s)) - X s £ (s))\ ds = 0(1 V T), and so ft H 5 (z 5 £ (s),0)ds = 0(1 V 
T). ' ' □ 

We now follow steps one through four from Section 14. H making modifications where 
necessary. 

Step 1: Reduction using Gronwall's Inequality. Now h s e {r/e) satisfies h s e {r/e) — 
/if (0) = e j T/£ H s {z s £ {s),0)ds. Define e*(r) =_e J q t/£ H s {z s e {s),0) - H s {h 5 £ {s))ds. It follows 
from Gronwall's Inequality and the fact that H'(-) G C 1 ({(5, h) : < 5 < Sq, h G V}) that 



sup \h s e (r/e) - h s (r)\ < 

Q<r<TAT* 



sup , 

, Q<t<TAT* 



r; ( r)| I e Lip ( 54 lv)^ = o | sup 

0<t<TAT« 



e 5 e (r) 



Step 2: A splitting according to particles. Next, 





H 5 {z,0) -H 5 (h) 



-k' 5 {X - Xl (z)) - 

o 



+ 



K ' 6 ( X2 ( Z ) - X) + 



-Wk' s (x 2 (z) - X) - 



and so, in order to show that sup 0<r<TAT i e*(r) = 0(e), it suffices to show that for i — 1,2, 



sup 

0<r<TAT/ 



sup 

o<t<tatS 



r/ 6 



8m^f i£ ( S ) 



., «')-"«W)l) + v; W , l4 (,),„ 



r/e 



8m^J e (s) 



T,(^( S ),^ £ ( S ),5) 



0(1). 



We only demonstrate that sup 0<T 



<TAT* 



<(.Y| (S) - „«(.))) + gj^^j. 



f/.S' 



(9(1). The other three terms are handled similarly. 
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Step 3: A sequence of times adapted for ergodization. Define the sequence of times 
t% £ inductively by t 5 0>e = wf{t > : tf^t) = 0}, t s k+1>e = mi{t > t% E : <p{Jf) = 0}. If e 
and 5 are sufficiently small and t 5 k+le < (T A T £ )/e, then it follows from Lemma 15.21 that 
1/uj < t 5 k+le — t k£ <uj. From Lemmas 15.21 and [5.41 it follows that 



sup 

0<r<TAT| 



<0(1) 



8mi^f e (s) 



; /e KiWW -. l( ,fw)) +FiWWi£Us)ii) 



ds 



fe+l,e — 



■k + l,e 



<(X s £ (s) - xi(4{ 8 ))) + 



8mi^f e (s) 



(19) 



Step 4: Control of individual terms by comparison with solutions along fibers. 

As before, it suffices to show that each term in the sum in Equation (|T9|) is no larger than 
0(e). Without loss of generality we will only examine the first term and suppose that 

4e = °> i- e - that </?L(o) = o. 



Lemma 5.5. Iftf >E < then sup < t < t « £ |zg(t) - z 5 £ {t)\ = 0(e). 

Proof. By Lemma l5~l h s (t) - h s £ (t) = /if (0) - h s £ (t) = -e f* H 5 (z 5 £ (s), 0)ds = 0(e(l V £)) for 
t > 0. 

Using what we know about the divergence of the slow variables, we find that 



¥>i,o(f) - < £ (i) 



1 



T^XUs), E s (s), S) T X {X»{8\ El(s), 5) 
0{e)ds = 0{e) 



+ 0{e)ds 



for < t < t\ £ . Lemmas 15.11 and 15.21 ensure the desired uniformity in the sizes of the orders 
of magnitudes. Showing that sup 0<<< ^ |<^2,o(0 ~~ V2,e(0| = ®( £ ) * s similar. 

From Lemma 1531 we find that t 1>£ = t lfi + 0(e) = Ti(X 5 , E%, 5) + O(e). Hence 



□ 



8miEf £ (s) 



, T^a), 
But when Xi(z 5 £ ) < X s £ - a, 



ds = 0{e) + 



8miM 



-ds = O(e) + J8miEl Q 



Kt(X»{8) - X X {4{8))) 



sig: 



n(v l (4(s)))K' s {X s £ (s)-x l (4(s))) 



and so 



k' 5 (X s e (s) - Xl (z 5 M))ds = -^rmElM - v /2m 1 ^ )E (tf i£ ) = 0{e) - ^/Sm^o- 
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Hence, 

,<f„ JsmtE'Js) 

as desired. 



6 Appendix to Section El 

Proof of Lemma 15. It 

Proof. For < 5 < 8/2, 



We only consider the claims about T\, and for convenience we take m\ = 2. Then 

r x ~ a ds r x/2 ds 



E 1 - U X { 8 ) J a JE X - K S {S) 

X/2-5 f s ds \ 2X-A5 t . f 1 ds 
l _ + 45 ' 



-E'l ia \jEi — k s (s) J y/Ei Jk- 1 {e 1 ) \/E\ — k(s) 

Define 

J k-^{e) \/ E - k(s) Jo VE-u 

Notice that (k -1 )'^) diverges as u — > + , while (E — u)" 1 ! 2 diverges as u — > but both 
functions are still integrable on [0,-E 1 ]. It follows that F(E) is well defined. Then it suffices 
to show that F : [£, k(0) - £] -> R is C 1 . 
Write 

«/2 -( K -i)>(u) , ^ -(k^'M 



= / "rj^-' du + / "r^-' du = F^E) + F 2 (£). 
Jo \l E - u Js/2 V E - u 

A standard application of the Dominated Convergence Theorem allows us to differentiate 
inside the integral and conclude that F\ G C°°([£, k(0) — £}), with 

FliE) = L 2(E-u)^ du - 
To examine F 2 , we make the substitution v = E — u to find that 

F 2 {E)= / '-)= J -dv. 



o 
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Using the fact that (« 1 ) / G C 1 ([S/2, k(0)]) and the Dominated Convergence Theorem, we 
find that F 2 is different iable, with 



- S/2 Jo 

Another application of the Dominated Convergence Theorem shows that F 2 is continuous, 
and so F 2 e C\[E,k(G) - £]). 
Thus 

IX r 

T X {X, E u S) = -= + 4S -E- 1 ' 2 + F^EJ + F 2 (E l ) 



/ E 1 

has the desired regularity. For future reference, we note that 

&Ti 2 &r x -X 

□ 

Corollary 6.1. For all 5 sufficiently small, the flow z^t) restricted to the invariant tori 
A4 C = {h = c} is ergodic (with respect to the invariant Lebesgue measure) for almost every 

ceU. 

Proof. The flow is ergodic whenever the periods T\ and T 2 are irrationally related. Fix 5 
sufficiently small such that = —XjE^ 2 + 0(5) < 0. Next, consider X, W, and E 2 fixed, 
so that T 2 is constant. Because T\ EC 1 , it follows that, as we let E\ vary, p- ^ Q for almost 
every E\. The result follows from Fubini's Theorem. □ 

Proof of Lemma 15.21 

Proof. For the duration of this proof, we consider the dynamics for a small, fixed value of 
5 > 0, which we generally suppress in our notation. For convenience, we take mi = 2. 

Let ip denote the map taking (X, W, x\, v%, x 2 , v 2 ) to (X, W, E\, E 2 , <pi, ip 2 ). We claim that 
ip is a C 1 change of coordinates on the domain of interest. Since E\ = vf + k$(xi) + k$(X — x\), 
Ei is a C 2 function of X\, V\, and X. A similar statement holds for E 2 . 

The angular coordinates (pi(xi, Vi, X) are defined by Equation ()14|1. We only consider (pi, 
as the statements for (p 2 are similar. Then (fi(xi, V\, X) is clearly C 1 whenever x\ ^ a, X — a. 
The apparent difficulties in regularity at the turning points are only a result of how the 
definition of tpi is presented in Equation (j!4jl. Recall that the angle variables are actually 
defined by integrating the elapsed time along orbits, and our previous definition expressed 
ip i in a manner which emphasized the dependence on x\. In fact, whenever |t>i| < VEi, 

j-l-J^( K fy(E 1 -v 2 )dv it Xl <5 

U + i Jo ( K s )( E ^ ~ v ) dv &xi> X -8. 

Here E\ is implicitly considered to be a function of x\,Vi, and X. One can verify that Dip 
is non-degenerate on the domain of interest, and so ip is indeed a C 1 change of coordinates. 
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Next observe that dtpi^/dt = 1/Ti, so Hadamard's Lemma implies that dipi jE /dt = 
1/Tj + 0(ef{8)). It remains to show that, in fact, we may take f(S) = 1. It is easy to verify 
this whenever x\ < X — 5 because dE\/dt = there. We only perform the more difficult 
verification when x± > X — 5. 

When x\ > X — 5, |t>i| < \f~E\ and E\ = v \ + ks{X — x\). From Equation (}2~Tj) we find 



that 



1 

^ + 
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2 ■ T 1{ X,E U S) • (O'W-O*- 
To find dipi/dt, we consider ipi as a function of vi,X, and so that 

difi difi dv\ dfi dX d(p% dE\ 
~~d7 = dv^~dt + ~dXltt + dE 1 dt ' 
Then, using Equations (J2H|) and (J22J), we compute 



(22) 



9^3 1 d^l 



<%! dt T\ v 5 

1/2 - <9Ti 



A{x- Xl ) 



dX dt 
dip i dEi 
dE x dt 



Tx dX 
1/2 - Vl 9T X 



eW) = eW 
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2 " T a ' 

1/2 -cp! 2 



+ 



Using that k' s (X - x x ) = kI^kT 1 ^ - t> 2 ))/5 = ^(k" 1 )'^ - v 2 ))" 1 , we find that 



9(^i 

'dE~ 1 ~dt 



eO 



1/2 - 
<5 



+ £0 



(OC^i - y 2 )cfo 



But here 1/2 — <p x is 0(5). See the proof of Lemma 15. Hl below. Thus the claims about dipi/dt 
will be proven, provided we can uniformly bound 



(^i - v 2 )dv. 



^'(Ex-vl 

Note that the apparent divergence of the integral as \v x 



> yE\ is entirely due to the 
fact that our expression for <p x from Equation ()22)1 requires \v x \ < \/E~\. If we make the 
substitution u = E\ — v 2 and let e = E\ — v 2 , then it suffices to show that 



sup sup 

£<E 1 <k(0)-£ 0<e<£i 



^)'(e] 



du 



< +oo. 



VEi-u 

The only difficulties occur when e is close to 0. Thus it suffices to show that 

1 



sup sup 

£<Ei<k(0)— £ 0<e<£/2 



is finite. But this is bounded by 
1 



sup 

0<e<£/2 



—IV 



(e) 



8/2 U^Yiu) 



-du 



(K-i)'(e) 



sup 

0<e<£/2 



8,2 (K- r )"M 



\JE\ — u 



du 



(S/2) 



which is finite because (k 
similarly. 



-ly/ 



-oo as e — > + . The claims about d^/dt can be proven 



□ 
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Proof of Lemma 15.31 

Proof. We continue in the notation of the proofs of Lemmas 15.11 and 15.21 above, and we 
set rrii — 2. Then from Equation (|22|) . we see that n' s (X — x{) = unless \(fi — 1/2| < 

ff I^i^YiEi ~ v 2 )dv = 5F(E 1 )/T 1 = 0{5). Dealing with (p 2 is similar. 

□ 
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